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fS| ■ We complete the program outlined in the paper of the author with A. Migdal and 

sum up exactly all the fluctuations around the instanton solution of the randomly large 
scale driven Burgers equation. We choose the force correlation function k to be exactly 
quadratic function of the coordinate difference. The resulting probability distribution 
satisfy the differential equation proposed by Polyakov without an anomaly term. The 
\q ' result shows that unless the anomaly term is indeed absent it must come from other 

ON ■ possible instanton solutions, and not from the fluctuations. 



In paper |1J we computed the probability distribution P(Su, r) of observing the velocity 
difference 5u at a distance r for the randomly large scale driven Burgers equation in the 
WKB approximation. The results obtained there coincided with the one obtained earlier 
in 0. Naturally, we were able to derive the probability distribution only for 5u ^> 0, the 
^ ; so-called "right tail" of P. We found 

P(Su,r) « exp I — — I (1) 

At the same time, the left tail, conjectured in [0 to be 

3 



P(5u, r) » r , - J(u 2 } < 5u < (2) 

\5u\ 2 

remained to be computed. 

To compute the left tail, we have to go beyond the leading instanton approximation. In 
this letter we show how to sum up all the fluctuations around the instanton found in H in an 
exact way. We find that the effect of the fluctuations around the instanton is equivalent to a 
simple quantum mechanics (found previously in 0) and therefore the summation is reduced 
to much simpler problem of solving that quantum mechanics. 



1 



As was discussed in [|1] the randomly driven Burgers equation is equivalent to the field 
theory with the action 

S = — i [ [i(ut + uu x — i>u xx ) H — f dxdydt fi(x, t)n{x — y)fi(y, t) (3) 



with coordinates u(x) and momenta fi(x), the probability distribution of observing the velocity 
profile u(x) being the wave function in the coordinate representation for (^). 

In [|T| the wave function in the momentum representation was constructed in the WKB, 
or instanton, approximation. Long range force-force correlation was considered, 

(f{x,t)f{y,t')) = 5{t-t')K{x- y ) (4) 

k is assumed to be a slow varying even function of x which behaves as 



k(x) « «(0) - —x\ |d<4/^=L (5) 
2 V 

and quickly turns into zero when |x| ^> L. 

While only special wave functions in the momentum representation were considered, let 
us for completeness give a more general instanton solution. It is rather easy to generalize the 
instanton found in |3J slightly to see that any correlation function (exp (/ dx\(x)u(x))) can 
be found as long as / dx\(x) = 0. The velocity retains its linear profile and the answer is 
justQ 

/ 2kq 



exp ^ dx\(x)u(x)j\ = exp 



dxxX(x) 



dxxX(x) 3> (6) 



The support of X(x) must lie within the ±L interval. 

The natural thing to do now would be to consider the fluctuations around the instanton 
solution, thus computing the corrections to the answer ([]). That was the program outlined 
in!. 

In principle, the next order correction to (@) is given by the determinant of the operator we 
derive by expanding the action (Q) around the instanton solution. In practice, to compute the 
determinant directly would be rather inconvenient. So we take a slightly different approach 
suggested in |TJ. Instead of computing the quantity^] 

Z[AoA(x)] = ^exp ^ dxXoX(x)u(x)^J (7) 

we introduce 

F = ^-\ogZ[X X(x)} (8) 

This quantity, as was shown in 0, is easy to expand around the instanton solution. We 
choose 

u = u inst + u (9) 
fi = AW + A* 



1 The construction of jp is equivalent to taking X(x) — Ao (6(x — Xq/2) — S(x + x /2)) 
2 We introduce an auxiliary parameter Aq which we can later set to 1. 
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u inst and /x inst being the instanton solution, and find 



/ dx\(x) f VuDJx u{x) exp [—Sj 

dx\{x)u inst (x) + — — — 7 ~a (10) 

J VuVfj, exp l—b J 

The action S is the expansion of the action (|3]), 

S = -i J p,(ut - uu xx ) + ^ J dxdydt p,(x,i)K(x - y)p,(y,t) + (11) 
d 



The last term of (|TTJ) is the interaction. We will refer to the action without the interaction 
term as S . 

The next step is to define the Green's functions of small fluctuations around the instanton. 
These are the Green's function of the free field theory provided for us by Sq. For example, if 
we define 

G(x, t; x', t') = i (/i(x, t)u(x' , t')) (12) 



o 



N{x,tyx',t') = (u(x,t)u(x',t')} 
M(x,t;x',t') = -(fl(x,t)fl(x' ',*')) 



o 



the brackets ( )o signifying the fact that we define these functions with the action So. We can 
easily find the equations these Green's functions satisfy, for instance 

N t+fa (umstN) - vN xx = - Jdy k(x- y) G{y, t; x', f) (13) 
G t + u inst ^ + vG xx = -i + 5(x - x')5{t - 

We could try to solve these equations, find the Green's functions and then start expanding 
(|l(J) in powers of the interaction. The perturbation theory (complete with Feynman diagrams) 
we get in this way is actually a perturbation theory in powers of 1/A. In that sense, it is much 
better defined than the perturbation theory we would get if we expanded the action around 
the zero value of the fields, as in the standard Wyld's technique (see 0). 

Using the perturbation theory around the instanton we could construct F as a function 
of Ao and then we could integrate it to reconstruct the quantity Z. 

However, to solve the equations ( |I3"D in the straightforward way is difficult. If, for instance, 
we used the expansion (||) we could try to expand the function iV in powers of x and x' as 
well, 

N(x, t; x', tf) = a(t, t') + b(t, t')xx' + c(t, t') {x 2 + x' 2 ) (14) 

and try to find the functions a, h and c. This approach would not work as a and c will 
diverge at t — > — oo. That shows that the function iV is more complicated than just a simple 
polynomial. 

That is why we are going to use a certain trick which will allow us to avoid solving (p^3|) . 
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We begin with writing down the instanton equations of motion again, 

u t + uu x - vu xx = -i J dy k(x - y)/J>(y) (15) 

fit + ufi x + v[i xx = (16) 

These equations, since they follow from the action (|3]), have at least three conservation 
laws. The first two are the generalized momentum and energy, the corresponding components 
of the energy-momentum tensor, 

E = i J dxfi(uu x — vu xx ) — — J dxdy ^l{x)k{x — y)[i{y) (17) 
dxjiu x 



They were originally found in |Q. In addition to these two, there is another conservation law 
following from the Galilean invariance of (^). The transformation laws are 

u(x) — > u(x — vt) + v (18) 
n(x) — > fi(x — vt) 

and the corresponding conserved quantity 

G = J dxfi-tP (19) 

If the generalized momenum P is zero (and it will be zero for all the solutions of ( pT5[ ) and 
(|16|) which fall to zero at t — > — oo), then even the quantity 

G = J dxn (20) 

is conserved and zero. The quantity G can be shown to be a canonical momentum conjugate 
to the constant component of the velocity field. It has an evident meaning of the center of 
mass of Burgers fluid (moving with the velocity given by the momentum P) . 

It is important to note that the Galilean symmetry is, in fact, spontaneously broken since 

u 2 \ = const (21) 

as was pointed out in . Spontaneous symmetry breaking is, as always, due to the boundary 
conditions in path integral. The natural boundary condition is u(t = — oo) = 0. If we take 
u(t = — oo) = v, we obtain another ground state of the system which can be taken to the first 
one by the symmetry transformation (TS). By the way, the instanton found in M satisfied 



the boundary condition u(t = — oo) = and therefore it was not invariant under (18); by 
applying (|T^) to it we could construct other instantons. We did not need to sum over all those 
instantons, though, because the symmetry is broken. Fortunately, this spontaneous symmetry 
breaking does not prevent us from exploring the consequences of the Galilean symmetry. 

We saw that the Galilean symmetry constraints the momentum conjugate to the constant 
component of the velocity which was the culprit of the unstability in the Green's functions 
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(compare with flUP). Could we find a way to constraint the integral over the momentum 
/ dxfi not only on classical solutions, but also in all the fluctuations? It turns out there is 
such a way 



Following H, we consider a slight modification of (18) allowing v to depend on time. This 



is no longer a symmetry of the action ([$]). Under the proposed transformation 

u(x) -> u(x-r(t)) + (22) 

(JjL 

fjb(x) — > fi (x — r(t)) 
dv 

r(0) = 0, — (t = -oo) = 

the action (^) changes as 

S ^ S -i J dt\^ J dx^ (23) 

Note that the transformation does not change the ground state due to the last condition in 

It turns out that there is a procedure first proposed in |3| (for a diferent reason) which is 
similar to the computation of the Faddeev- Popov determinant and which allows to integrate 
out the degree of freedom associated with r(t). We would like to give here the outline of this 
procedure. We begin with inserting into the functional integral / T>uT>fiexp(—S) the identity 

Vr5 (u(r) - J[u] = 1 (24) 

which serves as a definition of the quantity J[u], the Faddeev- Popov determinant. Then we 
shift the variables of integration according to ( [22] ) and integrate out r(t) to obtain the original 
path integral with the Faddeev-Popov determinant and two constraints. 



J VuVfiexp(-S) J[u) 6(u{0)) SU dx^j (25) 



To proceed further we need to compute the Jacobian J[u}. This Jacobian is just 

We have to understand that the operator whose Jacobian we are computing defines the time 
propagation backwards; r(0) = and then we find r(t) for t < 0. That allows us to compute 
the Jacobian using the standard methods (see, for example, 0) to get 



J[u] = exp 



9(0) / dt u x (0) 



(27) 



In this form the formula is rather ambiguous and depends on how we define the ^-function at 
0. The necessity to define the 9 function at is often encountered in path integral formalism. 



It is related to the way we discretize time in deriving (27) . By discretizing time we break 
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Galilean invariance and it is not evident that it will be restored in the limit when we take the 
discretization time to zero. We have to choose such discretization so as to preserve Galilean 
invariance in the continuous limit. 

Unfortunately, so far no one has managed to show by a direct computation what the correct 
choice for 6(0) should be. However, in the paper || an indirect procedure was developed 
which, when applied to the formalism in hand, supports the choice 6(0) = 1. Such a choice 
would be equivalent to taking the retarded discretization of Burgers equation and was already 
suggested in [|J. As we will see later, it is also consistent with the Polyakov's equation. 

Therefore, the Jacobian must be taken as 



J[u] = exp 



dt u x (0) 



(28) 



Now note that the two constraints precisely eliminate two unstable modes (compare with 
(p]))- The Jacobian is in fact an effective contribution of the modes we eliminated in 
this way. 

Now we can safely go on and find the Green's functions. The equations for the Green's 
function will be changed slightly to account for the constraints we have. For example, the 
equations ( pT3| ) will become (we use Mi nst = ax) 



d 

N t + a— (xN) - vN xx = 
dG 

Gt + ax— — h vG xx = —i 

ox ox 



-Jdy k(x - y) G(y, t; x\ t') + f(t; x\ t') 
d ^. N + [S(x - x') - 5(x)]5(t - f ) 



(29) 



the difference of the delta functions being necessary to ensure / dx G = and G(x' = 0) = 
and the function f(t; x', t') being chosen to ensure N(x = 0) = N(x' = 0) = 0. 

It turns out that unlike ( |T3| ) it is relatively easy to solve (p9[). We expand the function k 
in Taylor series and immediately see that the function iV is bilinear in x and x', 



N(x, t; x , t') = b(t, t')xx' 



(30) 



where the function b(t, t 1 ) coincides exactly with the correlation of fluctuations of the velocity 
strain. The functions a and c from fllil) disappear due to the constraints. Moreover, by a 
direct computation one can show that the perturbative expansion for the quantity Z defined 
in (0) coincides term by term with a perturbative expansion of the quantity 



\^exp 



crAo J dxxX(x) 



defined in a theory with the action 

Sqm = J dt \ip(a t + cr 2 ) ~ yP 2 + o\ 

This action follows from f ^5|) if we formally substitute u = a(t)x and p 
last term o comes out of the Jacobian (|28|). 



(31) 
(32) 

/ dx xfi(x). The 
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Sqm defines the motion in the potential^ 



a 4 



V{a) = — - 2a (33) 

This is the quantum mechanics discussed in [[| without the anomaly term (in the notations 
of 6=10)- 

This quantum mechanics was shown in || not to have a stationary state and gives an 
asymptotics to the (nonstationary) probability distribution 

P ( 6u ) ~ uTTT> Su < (34) 

This is certainly not what is observed in the numerical experiments, ||. And we note that 
the viscosity did not contribute to the calculations. 

This result is physically unsatisfactory. We would have expected to find the anomalous 
contribution to the Polyakov's equation due to nonzero viscosity. One possible explanation 
could be that the instanton found in ][]] is not the only possible instanton solution. If this is 
so, then the further program is to find the instanton solutions for Burgers problem which do 
not use the Taylor expansion of the force correlation function. Perhaps they will contain the 
shock waves combined with linear profiles and the fluctuations around them will produce 
us the anomaly term of and lead to the stationary state. 

This problem notwithstanding, we would like to conclude with saying that we produced an 
answer for the Burgers problem in the region where the shock waves are absent by summing 
up exactly the perturbations around the instanton found in 

The author is grateful to A. Migdal who was a collaborator in this project at its ear- 
lier stages, and to S. Boldyrev, A. Polyakov, E. Balkovsky, V. Lebedev, G. Falkovich and I. 
Kolokolov for numerous and important discussions which helped to find the correct interpre- 
tation to the results obtained in this paper. 
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